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It is shown that /-electron multipole is naturally defined as spin-charge one-electron density operator 
in the second-quantized form with the use of tensor operator on the analogy of multipole expansion of 
electromagnetic potential from charge distribution in electromagnetism. Due to this definition of multipole, 
it is possible to determine multipole state from a microscopic viewpoint on the basis of the standard linear 
response theory for multipole susceptibility. In order to discuss multipole properties of filled skutterudites, 
we analyze a seven-orbital impurity Anderson model by employing a numerical renormalization group 
method. We show our results on possible multipole states of filled skutterudite compounds. 

KEYWORDS: multipole, spin-charge density, filled skutterudites, numerical renormalization group 
method 



1. Introduction 

One of current research trends in material science is to 
synthesize new compounds which exhibit exotic phenomena 
concerning magnetism and superconductivity. Due to the im- 
provement of single crystal quality with the use of modern 
crystal growth techniques, essentially new phenomena have 
been also discovered even in known compounds. In particular, 
ordering of higher-rank multipole has been actively studied 
both from experimental and theoretical sides in the research 
field of strongly correlated /-electron systems. 1 ' 2 In general, 
owing to the strong spin-orbit coupling of / electrons, spin- 
orbital complex degree of freedom, i.e., multipole, is consid- 
ered to be active in /-electron compounds. However, when 
orbital degeneracy is lifted, for instance, due to the crystal 
structure with low symmetry, only spin degree of freedom of- 
ten remains. Namely, /-electron compounds crystallizing in 
the cubic structure with high symmetry are quite important 
for the research of multipole phenomena. 

In this sense, filled skutterudite compounds LnT4Xi2 with 
lanthanide Ln, transition metal atom T, and pnictogen X pro- 
vide us an ideal stage for the promotion of multipole physics, 
since this material group crystallizes in the cubic structure 
'■h . of Th point group. 3 Moreover, it is possible to synthesize 



many isostructural materials with different kinds of rare-earth 
and actinide ions, leading to the development of systematic 
research on multipole ordering. In fact, recent experiments 
in cooperation with phenomenological theory have revealed 
that multipole ordering frequently appears in filled skutteru- 
dites. For instance, a rich phase diagram of PrOs4Sbi2 with 
field-induced quadrupole order has been unveiled experimen- 
tally and theoretically. 4-6 Recently, antiferro retype higher 
multipole order 1 has been discussed for P1-RU4P12 7 ' 8 and 
PrFe4Pi2. 9-11 Except for Pr-based filled skutterudites, signs 
of multipole phenomena have been also found. In NdFe4Pi2, 
a significant role of quadrupole at low temperatures has been 
suggested from the measurement of elastic constant. 12 In 
SmRu4Pi2, a possibility of octupole order has been proposed 
from several kinds of experiments. 1317 

As mentioned above, theoretical research on multipole or- 
der has been developed mainly from a phenomenological 
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viewpoint on the basis of an LS coupling scheme for multi- 
/-electron state. Such investigations could explain several ex- 
perimental results and research activity along this direction 
will be still important in future. However, we strongly believe 
that it is also important to promote microscopic approach for 
understanding of multipole phenomena in parallel with phe- 
nomenological research. Based on this belief, the present au- 
thor and collaborators have made effort to develop a micro- 
scopic theory for multipole-related phenomena by exploiting 
a j-j coupling scheme. 18 In particular, octupole ordering in 
Np02 has been clarified by evaluating multipole interaction 
with the use of the standard perturbation method in terms of 
electron hopping. 19-21 We have also discussed possible multi- 
pole states of filled skutterudites by analyzing multipole sus- 
ceptibility of a multiorbital Anderson model based on the j-j 
coupling scheme. 22 "" 27 

However, it seems to be still difficult to understand intu- 
itively the physical meaning of multipole degree of freedom 
due to the mathematically complicated form of multipole op- 
erator defined by using total angular momentum. As men- 
tioned above, multipole is considered to be spin-orbital com- 
bined degree of freedom. In this sense, it seems to be natu- 
ral to regard multipole as anisotropic spin and/or charge den- 
sity. This point has been first emphasized in the visualization 
of octupole order in Np02- 19-21 The definition of multipole 
as spin-charge density has been briefly discussed, when we 
have attempted to clarify multipole state of heavy lanthanide 
filled skutterudite. 28 Owing to the definition of multipole in 
the form of one-electron spin-charge density operator, it has 
been possible to discuss unambiguously multipole state by 
evaluating multipole susceptibility even for heavy rare-earth 
compounds with large total angular momentum. 

In this paper, first we explain the definition of multipole 
as spin-charge density in the form of one-body operator from 
the viewpoint of multipole expansion of electromagnetic po- 
tential from charge distribution in electromagnetism. In or- 
der to determine the multipole state, we use the optimization 
of multipole susceptibility on the basis of the standard lin- 
ear response theory. To proceed with further discussion, in 
this paper we pick up an impurity Anderson model including 
seven / orbitals. We perform the calculation of multipole sus- 
ceptibility by using a numerical renormalization group tech- 



2 J. Phys. Soc. Jpn. 



Full Paper 



Takashi HOTTA 



nique. Then, we review our recent results of heavy rare-earth 
filled skutterudites for Ln=Gd^Yb. We also show our new 
result for multipole susceptibility of Sm-based filled skutteru- 
dites with some comments on the effect of rattling. Finally, 
we summarize this paper and briefly discuss a future problem 
of multipole theory. Throughout this paper, we use such units 
as h=ks=l and the energy unit is set as eV. 

2. Multipole as Spin-Charge Density 

In general, multipole is a concept to express the degree of 
deviation from spherical symmetric structure. A well-known 
example can be found in the multipole expansion of electro- 
magnetic potential from charge distribution. The potential is 
formally expanded by the spherical harmonics Ylm and both 
electric and magnetic multipole moments are defined from the 
expansion coefficients. 

Here we emphasize that the multipole expansion is also 
applicable for the consideration of electromagnetic field pro- 
duced by electrons. 29 For instance, Kubo and Kuramoto have 
discussed the multipole expansion of the vector potential from 
local electrons, in order to estimate the internal magnetic field 
from octupole moment in Ce x Lai_ x B6. 30 

Based on the above background, now we consider the def- 
inition of multipole as /-electron density operator X, which 
is generally expressed in the second-quantized form as 



X ^ ' X ma m i a i f ma fm'a' i 



(1) 



ma.ra' er' 



where f ma is the annihilation operator for / electron with 
spin a and z-component m of angular momentum £=3. 
Throughout this paper, we define a=l (—1) for up (down) 
spin. The form of the coefficient X mt7iTn ' a i is related to the 
definition of multipole. As easily understood, X denotes total 
charge, i.e., monopole, for 

-^monopolc c c 



(2) 



Intuitively, we can understand that /-electron multipole de- 
notes the deviation from the total charge operator. In the stan- 
dard definition, dipole denotes total angular momentum J. 
Thus, X for dipole is given by 

-v^dipolc _ tq _ c rot , not c /o\ 

^ma,m'<j' ~ J ma-,m'cr' ~ °<rv' ^mm 1 ' °aa' °mm' , (->) 

where a indicates Cartesian component, L denotes angular 
momentum operator for £=3 and S indicates spin operator. 
Then, how do we define higher-order multipoles? 

In the multipole expansion of potential in electromag- 
netism, higher electric and magnetic multipole moments ap- 
pear in the coefficients of the expansion by the spherical har- 
monics Ylm with larger angular momentum. In group theory, 
Ylm is defined by the basis of irreducible representation 
of the rotation group R, expressed as 



RYlm = J2Y LM ,D 



(L) 

MM'- 



(4) 



M> 



On the analogy of the multipole expansion, for /-electron 
multipole operator, we exploit a concept of spherical tensor 
operator in the quantum mechanics of angular momentum. 31 
When we consider the rotation of operator T, we obtain a 
set of operators f W = {f 9 (fe) } with (2k + l)-components 



(q = —k, — k + 1, • • • ,k—l,k), given by 



(5) 



q' 



* (k) 

Namely, Tq is transformed like a basis of irreducible repre- 
sentation for the rotation. Such fq k ^ is called spherical 
tensor operator of rank k. 

" (k) 

The spherical tensor operator T g for / electron is ex- 
pressed in the second-quantized form as 

Erp(k.q) j-f p 
1 mo.m'a' J ma I'm'"' J 



f(k) 

q 



(6) 



mcr,m' a' 
-,(k,q) 



where the coefficient Tl,„'., is calculated as follows. First 
it is convenient to change the /-electron basis from (m, a) to 
(j, /i), where j is the total angular momentum and [i is the 
z component of j. Note that j takes 7/2 and 5/2 for / elec- 
trons. For a certain value of angular momentum j and its z- 
component /i, the matrix element of spherical tensor operator 
is easily calculated by the Wigner-Eckart theorem as 

0>|7f W) - {j ^f (j»\j»'kq), (7) 

where (JM\J'M'J"M") denotes the Clebsch-Gordan co- 
efficient and (,;?||T( fc )||j) is the reduced matrix element for 
spherical tensor operator, given by 



U\\T {k) \\j) 



1 / (2j + fc + 1)! 

2k" 



(2j - fc)! 



(8) 



Note that fc < 2j and the highest rank is 2j. The coefficient 
v , is obtained by returning to the basis of (m, a) from 



n(k,q) 



~ ma,m'<7- 



(j, n). The result is given by 

-n(k,q) 



T^,, = V (Mo*,/*,) 
x (j^\lms^)(jij,'\£m's^-), 



(9) 



where £=3, s=l/2, j=£±s, and fj, runs between — j and j. 

We note that it is not necessary to take double summations 
concerning j in eq. (9), since the matrix representation of to- 
tal angular momentum J is block-diagonalized in the (j, y)- 
basis. We have checked that the same results as eq. (9) are 
obtained when we calculate higher-order multipole operators 
by following the symmetrized expression of multiple products 

QfJ 31,33,34 

Thus far, we have implicitly assumed /-electron density 
in an isolated ion, but in actuality, rare-earth ions are in the 
cubic crystal structure. Then, it is convenient to change from 
spherical to cubic tensor operators, given by 



(10) 



where fc is a rank of multipole, an integer q runs between — fc 
and fc, 7 is a label to express Oh irreducible representation, 

f k) 

and G\/ q is the transformation matrix between spherical and 
cubic harmonics. Throughout this paper, we use the cubic ten- 
sor operator as multipole. 

It should be noted here that multipoles belonging to the 
same symmetry are mixed in general, even if the rank is dif- 
ferent. In addition, multipoles are also mixed due to the effect 
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of crystalline electric field (CEF) of Th point group. Namely, 
the /-electron spin-charge density should be given by the ap- 
propriate superposition of multipoles, expressed as 



X 



(ii) 



In order to determine the coefficient p\ , we evaluate the mul- 
tipole susceptibility in the linear response theory. 24 Namely, 

(k) 

p\ is determined by the normalized eigenstate of suscepti- 
bility matrix, defined as 

e -Ei/T _ e -E } /T 



1 



Xk-y,k'-/' 



Ej — Ei 



(12) 



where Ei is the eigenenergy for the i-th eigenstate |i) of the 

(k) 

Hamiltonian H of the system, T is a temperature, p\ — 
J2i e~ Ei t T (i\T~^ \i) /Z, and Z is the partition function given 
byZ = E i e-^ T . 

When we express the multipole moment as eq. (11), we 
normalize each multipole operator so as to satisfy the or- 

thonormal condition Tr{T^ fe) f^' ) }= 4fc'<5 77 '- 35 We note that 
the multipole susceptibility is given by the eigenvalue of 
the susceptibility matrix. Note also that the susceptibility for 
4u multipole moment does not mean magnetic susceptibil- 
ity, which is evaluated by the response of magnetic moment 
i+2S,i.e., J+S. 22 ' 23 

3. Model and Parameters 

In the previous section, we have defined multipole as spin- 
charge density in the form of one-body operator. In order to 
determine the multipole state, we have explained a method 
to optimize the multipole susceptibility based on the linear 
response theory. To proceed with further discussion, it is nec- 
essary to set the Hamiltonian H. For /-electron systems, it 
is desirable to treat a seven-orbital periodic Anderson model, 
but at least at present, it seems to be a heavy task to analyze 
such a multiorbital model. To reduce the task in calculations 
by keeping essential physics, we can take two ways. 

One is to decrease the number of relevant / orbitals in the 
periodic system. As for the research along this line, here we 
briefly introduce the theoretical study on octupole order in 
Np0 2 . 19-21 After the discussion about the CEF states of ac- 
tinide dioxides, two relevant Tg orbitals have been extracted. 
Then, on the basis of a j-j coupling scheme, the r 8 orbital de- 
generate Hubbard model has been set as an effective Hamil- 
tonian for Np02- With the use of the standard perturbation 
theory in terms of /-electron hopping, effective interactions 
between multipoles in adjacent sites have been evaluated on 
an fee lattice. Due to the combination of exact diagonalization 
and mean-field theory, it has been concluded that the ground 
state has longitudinal triple-q 5u octupole order, consistent 
with experimental facts. 

Another way is to keep all seven / orbitals, but to con- 
sider an impurity Anderson model. Even if the number of lo- 
cal degree of freedom is increased, it is possible to solve the 
impurity Anderson model with the use of a numerical renor- 
malization group method. In this paper, we adopt this way to 
discuss multipole state of filled skutterudites. 



The seven-orbital Anderson model for filled skutterudites 
is given by 

H = J2 £kC L Ck ^+ {Vmcijma +h.C.) + Hl oc , (13) 

where is conduction electron dispersion and c^a is the an- 
nihilation operator for conduction electron with momentum 
k and spin a. The second term indicates the hybridization 
between conduction and / electrons. For filled skutterudites, 
the main conduction band is given by a u , constructed from 
p-orbitals of pnictogen. 36 Note that the hybridization occurs 
between the states with the same symmetry. Since the a u con- 
duction band has xyz symmetry, we set V%=— V-2=V and zero 
for other m. We fix V as V=0.05 eV and a half of the band- 
width of a u conduction band is set as 1 eV. 

The third term in H is the local /-electron part, given by 



H\ oc = H so + Hi n t + He 



EF- 



(14) 



Note that the chemical potential is appropriately changed to 
adjust the local /-electron number n. The spin-orbit coupling 
term H so is given by 



H so — A 



£ 

m.a.m'a' 



m.a.m 



(15) 



where A is the spin-orbit coupling and f ma denotes the anni- 
hilation operator for / electron with spin a and angular mo- 
mentum m(=— 3,- • - ,3) and a=+l (—1) for up (down) spin. 
The matrix elements are expressed by Cm,±i,m,±i=i™/2, 
Cm±i,=Fi,m,±i=\/l2 — m(m ± l)/2, and zero for the other 
cases. The second term if; nt indicates the Coulomb interac- 
tions among / electrons, expressed by 

Hi n t~ ^ ' E ni ^ m2 , m:i , m4 f^ nirTi f^ n2a . 2 f m3r j 2 f mi( j 1 , (16) 



mi 
CT 1 > CT 2 



where the Coulomb integral / is expressed by the combination 
of Slater-Condon parameters, F°, F 2 , F 4 , and F 6 . 37 
Finally, the CEF term Hcef is given by 



-HCEF - Sm,m'/m<7/n 



(17) 



m.ra' .a 



where B m . m ' is determined from the table of Hutchings for 
angular momentum J=f=3, 38 since we are now considering 
the potential for / electron. For filled skutterudites with T h 
symmetry, 3 B m m i is expressed by using three CEF parame- 
ters B\, £?g, and B\? % Following the traditional notation, we 
define 



B 2 



Wx/F(4), 

W(l-\x\)/F(6), 

Wy/F\6), 



(18) 



where x and y specify the CEF scheme for Th point group, 3 
while W determines an energy scale for the CEF potential. 
Concerning F(4), F(6), and F*(6), we choose F(4)=15, 
F(6)=180, and F*(6)=24 for^=3. 38 - 39 

For further discussions, it is necessary to set the parame- 
ters. First let us consider the Slater-Condon parameters F k 
(k=Q, 2, 4, and 6) and spin-orbit coupling A. We will discuss 
CEF parameters later. In order to determine F k and A, we 
try to reproduce the / 2 excitation spectrum of Pr 3+ ion. 40 ' 41 
Here we consider twelve excited states, which are labeled as 
3 H 5 , 3 H 6 , 3 F 2 , 3 F 3 , 3 F 4 , 4 G 4 , 1 ^ 2 , 3 P , 3 Pi, 3 P2, *h, and 



4 J. Phys. Soc 

0.015 
0.010 



». 0.005 
>, 

£ 0.000 
c 

CD 



Jpn. 



Full Paper 



Takashi HOTTA 



-0.005 



-0.010 



-0.015 
-1 

0.004 



(a) W=-4x10" 4 | 


1 1 


y=0.3 


- 


n=2 ! 




\- 23 | 




r +(D 1 

^ r" 






/ r +(2) 


: — - — -~^jy'> 









-0.5 



0.0 

x 



^ 0.002 



£ 0.000 



LU 

o 



-0.002 



-0.004 



-1.0 



-0.5 



0.0 

x 



0.5 



1.0 



0.5 



1.0 



0.020 



> 0.010 

>, 

c 0.000 

CD 



LU 

o 



-0.010 



-0.020 



(a) W=-4x10" 
y=0.3 
n=7 



octet (J=7/2) 



-1.0 



-0.5 



1 1 1 1 

(b) W=-4x10" 4 


- ■ 1 1 




y=0.3 






n=5 




> 






05, 






>. 












CD 






C 






CD 






LL 






LU 


■ i 


.i i 


o 



0.020 



0.0 

x 



-0.010 



■0.020 



(b) W=-4x10 
y=0.3 
n=10 



0.5 



1.0 




Fig. 1. CEF energy levels for (a) n=2 and (b) n=5. 



Fig. 2. CEF energy levels for (a) n=l and (b) n=10. 



1 5o- By diagonalizing H so + H lnt , we obtain eigenenergies 
Ej and define excitation energy SEj=Ej — Eq, where Eq is 
the ground state energy for 3 H^. We determine F k and A so as 
to minimize the sum of the square of energy difference, given 

b y A = T,]ii( SE T p - 5e j) 2 - where 5E T P indicates ^ 

excitation energy experimentally determined. Note that F° 
cannot be determined in the present simple procedure, since 
it appears only as the offset of the energy level when local 
/-electron number is assumed to be unchanged. In order to 
determine the value of F° itself, it is necessary to resort to the 
first-principle calculation, but it is out of the scope of this pa- 
per. Here we simply set by hand F°=10 eV as a typical value 
for rare-earth ion. 

After the minimization of A, we obtain F 2 =8.75 eV, 
J F 4 =6.60 eV, F 6 =4.44 eV, and A=0.095 eV. These values are 
expected to be reasonable, since the experimental value of A 
is 0.094 eV. For other lanthanides, we use the same values 
for the Slater-Condon parameters, since there is no reason to 
change two-body interactions, even when the electron number 
is changed. On the other hand, since the spin-orbit interaction 
is considered to be sensitive to ion radius, we use experimen- 
tal values of A for each lanthanide such as A=0.144 eV (Sm), 
0.180 eV (Gd), 0.212 eV (Tb), 0.240 eV (Dy), 0.265 eV (Ho), 
0.295 eV (Er), 0.326 eV (Tm), and 0.356 eV (Yb). 42 

Concerning CEF parameters, we emphasize that electro- 
static potentials act on one /-electron state. In principle, it 
is not necessary to change the CEF parameters when we sub- 
stitute rare-earth ion, as long as we consider the same crys- 
tal structure. Namely, it is enough to determine the CEF pa- 
rameters for some filled skutterudite compound. First we set 
W=— 0.4 meV and y=0.3, which are considered to be typical 



values for filled skutterudites. However, the CEF state is dras- 
tically changed by x. From the CEF energy level for n=2, as 
shown in Fig. 1(a), we choose £=0.3 so as to reproduce quasi- 
quartet CEF scheme of PrOs4Sbi2. 43 ^ 5 We use such CEF pa- 
rameters for other rare-earth ions. Since we discuss multipole 
states for Sm-, Gd-, and Ho-based filled skutterudites, here 
we show the CEF energy schemes for n=5, 7, and 10. Results 
for other cases are found in Ref. 23, although we used dif- 
ferent parameters there. Namely, the CEF energy schemes are 
not changed drastically, as long as we use realistic values for 
Coulomb interaction and spin-orbit coupling. 

In Fig. 1(b), we show the CEF energy levels for n=5. There 
appear doublet and T7 7 quartet states for x < and 
x > 0, respectively. At £=0.3, we obtain T^ 7 quartet ground 
state. From the specific heat measurements for S111RU4P12 and 
S111OS4P12, the CEF ground state has been concluded to be 
Tgy quartet. 46 Recently, magnetization measurement has been 
performed and the observed anisotropy has confirmed the Tq 7 
quartet ground state in SmOs4Sbi2. 47 On the other hand, the 
CEF energy scheme can be changed even among the same 
Sm-based filled skutterudites. In fact, for SmFe4Pi2, dou- 
blet ground state has been suggested. 46 ' 48 We have explained 
that such conversion of the CEF ground state at n=5 occurs 
due to the balance between Coulomb interaction and spin- 
orbit coupling. 27 

In Figs. 2(a) and 2(b), we show the results for n=7 and 10, 
respectively. For the case of Gd 3+ with n=l, the CEF ground 
state is well described by L=0 and J=S=7/2, but it is almost 
independent of x, since the CEF potentials for L=Q provide 
only the energy shift. Note, however, that the CEF state is 
given by the mixture of the LS and j-j coupling schemes. 
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This point will be discussed later again. For the case of Ho 3+ 
with n=10, a remarkable point is that r^ 3 doublet and 
triplet states are almost degenerate in the wide range of the 
values of x. In particular, around at £=0.3, the ground state is 
easily converted. Thus, it is necessary to discuss carefully the 
case of n=10. 

4. Numerical Results 

In order to evaluate the susceptibility matrix, we employ a 
numerical renormalization group (NRG) method, 49 in which 
momentum space is logarithmically discretized to include ef- 
ficiently the conduction electrons near the Fermi energy. In 
actual calculations, we introduce a cut-off A for the logarith- 
mic discretization of the conduction band. Due to the limita- 
tion of computer resources, we keep M low-energy states. In 
this paper, we set A=5 and Af=4500. Note that the temper- 
ature T is defined as T=A~ ( - N ~ 1 ^ 2 in the NRG calculation, 
where N is the number of the renormalization step. 

First let us summarize our recent results for multipole state 
of heavy lanthanide filled skutterudites. 28 For Ln=Gd, we 
have found the effect of quadrupole moments due to the de- 
viation from the LS coupling scheme. For Ln=Ho, when the 
CEF ground state is r^ 3 doublet, the exotic state dominated 
by 2u octupole moment has been observed. For Ln=Tb and 
Tm, the CEF ground state is Tf singlet and we have found no 
significant multipole moments at low temperatures, although 
we cannot exclude a possibility of antiferro lg+2g ordering. 
For Ln=Dy, Er, and Yb, the CEF ground state is doublet, 
and the dominant moment is the mixture of 4u and 5u. Among 
these cases, we pick up the results of Ln=Gd and Ho, and re- 
view interesting possibilities for the multipole states. 

For the case of Gd 3+ ion, the dominant multipole com- 
ponent is 4u dipole and the secondary components are 3g 
and 5g quadrupoles. Note that after the partial screening 
due to conduction electrons at extremely low temperatures, 
quadrupole moments disappear, but it seems to observe sig- 
nificant contribution from quadrupole moments. In the LS 
coupling scheme, the / 7 state is specified by J=S=7/2 and 
L=Q. Namely, at the first glance, we do not expect the ap- 
pearance of quadrupole moment. However, we should note 
that actual situation is always between the LS and j-j cou- 
pling schemes. Namely, some finite contribution of the j-j 
coupling scheme is included in the ground state. If very large 
A is assumed, first j=5/2 sextet is fully occupied and then, 
one / electron is accommodated in j=7/2 octet. Thus, this 
state can be multipole-active. Recently, it has been observed 
in GdRu4Pi2 that 101 Ru NQR frequency exhibits temperature 
dependence below a Neel temperature Tn=22K. 50 This may 
be interpreted as the effect of quadrupole due to the deviation 
from the LS coupling scheme. 

For the case of Ho 3+ ion, the CEF ground state is triplet 
at £=0.3, but the first excited state is r^ 3 doublet with very 
small excitation energy such as 10~ 5 eV. Thus, as shown in 
Fig. 2(b), Ho-based filled skutterudite is considered to be in 
the quasi-quintet situation. 23,51 For x=0.3, we have found sev- 
eral kinds of multipoles, but the dominant one is always given 
by the mixture of 4u and 5u from dipole, octupole, dotriacon- 
tapole, and octacosahectapole. 

Here we note that the CEF ground state is fragile for n=10. 
If we slightly decrease x, the ground state is easily changed. 
Then, we have evaluated multipole susceptibility for £=0.25 




Fig. 3. (Color online) Multipole susceptibility for n=5. 

with r^j doublet ground state. The multipole states at high 
temperatures are similar to those for £=0.3, but at low tem- 
peratures, we have found 2u multipole state, expressed as 

(3) (7) 

p 2u =0.955 and p 2u =0.297. Namely, the main component is 
2u octupole, but there is contribution of 2u octacosahectapole. 
From the elastic constant measurement for HoFe4Pi2, some 
anomalous features have been discussed. 51 It seems interest- 
ing to reexamine experimental results from the viewpoint of 
2u octupole state. 

Now we show our new result for the case of n=5, corre- 
sponding to Sm-based filled skutterudites. As mentioned in 
the introduction part, SmRu4Pi2 has attracted much attention 
from a possibility of octupole ordering. It is meaningful to 
discuss possible multipole state in the present calculation. As 
shown in Fig. 3, we find that 3g quadrupole is dominant with 
p 3 ^=0.995 and p 3 ^=0.097. The secondary component is ex- 
pressed as 4u+5u, but the contribution from 5u is found to 
be less than one percent. Note that 4u and 5u are mixed due 
to the Th symmetry, although such mixing can be included 
neither in the LS nor j-j coupling scheme. It is an advan- 
tage of the present calculation to take into account the 4u-5u 
mixing correctly. The third and fourth components are, re- 
spectively, expressed by 4g+5g and 2u. The fifth component 
is also expressed by 4u+5u, which includes significant contri- 
bution from 5u octupole. 

It is difficult to conclude the ordering of multipole only 
from the present results, but the local multipole moment with 
significant weight is considered to be a candidate which is or- 
dered in actual system. Then, if we discuss relevant octupole 
moment within the present result of the electronic model, 4u 
octupole seems to be the best candidate, since the secondary 
component includes 4u octupole moment and the contribution 
of 5u octupole is very small. 

Throughout this paper, we have not included the effect of 
phonons, but in the filled skutterudite structure, anharmonic 
local phonon, i.e., rattling, has been considered to play some 
roles to determine electronic properties. When we have con- 
sidered that the dominant phonon mode in filled skutterudites 
is Jahn-Teller type with E g symmetry, 52 the multipole state 
for the case of n=5 has been found to be characterized by 
the mixture of 4u magnetic and 5u octupole moments with 
clear difference between longitudinal and transverse modes. 53 
If such symmetry lowering can be detected in experiments, 
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we should consider seriously the possibility of 5u octupole. 
We believe that it is an important test of 5u octupole ordering 
scenario for the Sm-based filled skutterudite. 

5. Summary and Comment 

In summary, we have defined the multipole as spin-charge 
one-electron density operator by using the cubic tensor op- 
erator. Based on the linear response theory, we have deter- 
mined the multipole state so as to maximize the multipole sus- 
ceptibility. Then, we have discussed possible multipole state 
of filled skutterudites by evaluating the multipole suscepti- 
bility of the seven orbital Anderson model with the use of 
NRG method. After the review of the results for Ln=Gd~Yb, 
we have shown the result for Sm-based filled skutterudites. 
Within the impurity Anderson model, quadrupole moment has 
been found to be dominant. We have observed that the sec- 
ondary component is magnetic with significant contribution 
of 4u octupole. Thus, when we simply ignore the effect of rat- 
tling, 4u octupole seems to be relevant in the Sm-based filled 
skutterudite. It is consistent with recent result on the specific 
heat measurement of SmRu 4 P 12 . 17 

Finally, we briefly comment on the multipole expansion of 
/-electron spin-charge density. In this paper, the spin-charge 
density has been expressed by the linear combination of cu- 
bic tensor operators in the cubic crystal structure. Unfortu- 
nately, it does not provide the general multipole expansion of 
/-electron spin-charge density, since we are restricted only in 
the cubic symmetry. It is one of future issues to complete the 
multipole expansion of /-electron spin-charge density. 
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